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I. INTRODUCTION 



The recent decade has witnessed extensive research 
efforts directed at localized surface plasmons (LSP's), 
i.e resonant charge-density oscillations confined to sub- 
wavelength metal structures, such as particles^, holes^, 
shells^ or rods^. LSP's are important because they lead 
to strongly enhanced optical absorption and scattering 
cross-sections, and because they readily couple to op- 
tical far-fields, unlike the ordinary surface plasmons of 
extended metal surfaces. In the small particle regime, 
the LSP resonance wavelengths of a single nanostruc- 
ture are uniquely determined by its shape and dielectric 
function, and by the optical constants of the embedding 
medium. Thus, the color of the nanostructure can be 
tuned over an extended wavelength-range, including most 
of the visible and infrared regions in the case of silver 
or gold structures, through a variation in shape. More- 
over, the color can be used to "sense" the optical proper- 
ties of the surrounding. In addition, excitation of LSP's 
leads to strongly enhanced and localized fields in certain 
regions near the metal surface, and these fields can be 
used to amplify various molecular cross-sections. A bi- 
sphere system, for example, support resonances for which 
the field is concentrated to the gap between the spheres, 
which is crucial for single-molecule surface-enhanced Ra- 
man scattering (SERS)^, and a sharp point or protru- 
sion support LSP-enhanced fields at its apex, which can 
be utilized for various types of near-field microscopy^. 
Sensing applications rely on the fact that it is only the 
optical properties of the material within the zone of high 
field-enhancement that strongly affect the LSP resonance 



condition/ The LSP spectrum, which can be measured 
through far-field extinction or Rayleigh scattering spec- 
troscopy, can therefore be used for ultrasensitive sensing 
applications aimed, for example, at quantifying various 
biomolecular recognition reactions. 

As is known from a range of experiments and cal- 
culations on nanostructured or flat metal surfaces cov- 
ered with chromophores, such as dye molecules, the in- 
teraction between a plasmonic structure and its dielec- 
tric environment can be very strong if the dielectric 
possesses an optical resonance that is degenerate with 
a surface plasmon resonance jnode^S^iiiiSiiiiiiSiiSiii 
This strong coupling also forms the basis for impor- 
tant applications in molecular spectroscopy, in particular 
surface-enhanced resonance Raman scattering (SERRS) 
and surface-enhanced fluorescence. ^^'^^'^^ For atoms and 
molecules confined in microcavities similar strong cou- 
pling can occur In many applications, it is advanta- 
geous to bind the molecular layer to the metal via specific 
functional groups, for example thiols, in which case the 
molecular transition dipole moment will have a more or 
less well-defined orientation relative to the local fields 
generated at plasmon resonance. This will in turn affect 
the degree of coupling between a molecular resonance and 
a plasmon. 

In this work, we present an analytical theory for calcu- 
lating the optical properties of a sub-wavelength metal- 
lic ellipsoid, the prototypical example of a nanostruc- 
ture supporting tunable LSP's, covered by an optically 
anisotropic molecular layer. We then utilize this for- 
malism for investigating strong coupling effects in nano- 
plasmonics. The quasi-static theory of dipolar plasmons 
for ellipsoids has played an important role in nano-optics. 
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as it can be used to model a wide range of nanoparti- 
cle shapes of practical interest, including rods, spheres, 
oblate and prolate particles, accessing a broad range of 
frequencies. The extension of this theory to ellipsoidal 
cores with anisotropic coatings means that it is now pos- 
sible to also investigate the effect of molecular orientation 
relative to the core surface. 



II. THEORY 

While the polarizability of a metallic particle with and 
without coating has been studied in the past, and can 
be found in text books, to our knowledge, the effect of 
anisotropy, ie., molecular orientation, has not been an- 
alyzed earlier in the context of plasmonics. Here, we 
recapitulate recent analytic results^^'^^ for the dipolar 
polarizability of an ellipsoid with an anisotropic coating 
(the coating dielectric function being different parallel 
and perpendicular to the coating normal, see Fig. ^), 
and combine these results with realistic microscopic di- 
electric functions for the metallic nanoparticle and the 
coating. 

The system we have in mind is depicted in Fig. 
We consider a coated ellipsoidal particle in an exter- 
nal electric field Eq^ with field-component Eoy in the 
'^-direction {v = x^y^z). The frequency of the exter- 
nal field is uj = 27rc/A, where c is the speed of light and 
A is the wavelength. The principal semiaxes of the inner 
and outer ellipsoids are and by with ax > ay > az 
and bx > by > bz. The shape of the coated particle is 
completely specified by the ellipticity of the inner sur- 



face 
face 



a^/a^, the ellipticity of the outer sur- 



for a sphere) and 



s = (a^ — ay)/{a'l — a^); ef , and s are all in the range 
[0,1]. The coating thickness is determined by the relative 
coating volume Ay = Koat/'K), where Vq = 4:7rbxbybz/3 
is the total ellipsoidal volume and Koat = Vq — Vi is the 
coating volume {Vi = Airaxayaz/S is the volume of the 
inner ellipsoid). For the case of a thin coating the ellip- 
ticities are related: e? = e^(l -hJ), where S is the relative 
coating thickness parameter (see Refs. I23L 24): the rela- 
tive coating volume (a dimensionless entity) introduced 
above can be written in terms of the parameter S accord- 
ing to Ay = Fo6/ 2, where Fo = l^{l-sel)-^^{l-el)-^ 
(see Ref. The "material" properties of the coated 

nanoparticle enters through the relevant dielectric func- 
tions: We denote the dielectric function of the inner el- 
lipsoid by £in{^)' The coating has dielectric function 
ec\\{uj) in the normal direction and ec±(^) in the tangen- 
tial direction. The dielectric function of the surrounding 
medium is assumed to be real and frequency independent 
and is denoted by Sm- The entities above completely de- 
termine the electromagnetic response of a coated ellip- 
soidal particle, i.e., determine, for instance, the particle 
polarizability a^y and the electric field distribution in 
and around the coated ellipsoid j^^*^ 

The full expression for the dipolar polarizability ayy 
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FIG. 1: a) A cut through a coated ellipsoidal nanoparti- 
cle. The principal semiaxes perpendicular to the paper are ay 
and by respectively. The relevant dielectric functions ein(cj), 
£c||(^), £c±{(^) and Sm are also indicated, b) Extinction spec- 
trum, where aw = v = x,y,z. The dash-dotted 
curves are for an uncoated metallic ellipsoid, where the red 
(blue) curve is for the external field along the x-axis (z-axis). 
The red solid (x-direction) and blue solid (z- direct ion) curves 
correspond to the case of an ellipsoidal molecular layer {with- 
out an interior metallic region, i.e. we set = 1) with a 
resonance parallel to the layer normal. Here, this situation 
is modeled by taking ^q\\/vq — 0.2, 7oo||/'^o = 0.02, and 
^q^/vq — ^oo^/vq — 0.01 (corresponding to effective oscilla- 
tor strength Co|| ~ 2.26, and large frequency refractive index 
n|| ^ 1.06, see text). The red dashed (x-direction) and blue 
dashed (z-direction) curves are for the case that the molecular 
layer has a resonance perpendicular to the layer normal. We 
take 7o±/'i^o 0.1 and ^oo±/vq 1oo\\/vq 0.01 
(so that co± ~ 1.13 and n^ ^ 1.07). The following parame- 
ters were used: el — 0-8, s = 1, Ay — 0.08, Tp/ujp — 0.05, 
^o||/^P — ^Q±l^p — 0.5, T\\/ujp — V^/ujp — 0.03 and Sm — 1. 
The extinction cross-sections are shown in units of the total 
volume Vo- The solid and dashed vertical bars indicate the 
plasmon and molecular resonance frequencies, respectively. 
Notice the small peak heights of the molecular layer reso- 
nances (here enlarged by a factor 50) compared with the plas- 
mon resonance heights. 



(defined through py = 4:7reoem<^yyEoy^ where py is the 
^-component, v = x^y or of the induced dipole mo- 
ment and eo is the permittivity of free space) in terms 
of the geometric and dielectric entities above is given in 
Refs. l23l and 24: A knowledge of ayy requires a function 
Hy{u,s;q)^ satisfying Heun's equation2^^, evaluated at 
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= e^, surfaces. The 



the inner, q = ef^ and the outer, q 
shape parameter s was introduced above. The anisotropy 
in the dielectric function of the coating {ec\\ and ec±) en- 
ters through u defined as 



u± 



(1) 



where = £cfi{^)/£m =^ ID- ^^e case of 

spheroids (two of the principal axes are equal) Hy{u, s; q) 
is expressible in terms of hypergeometric functions, which 
are available in standard mathematical numerical pack- 
ages, such as MatLab or Mathematica. For general el- 
lipsoids Hyiu^ s] q) is straightforwardly generated using a 
recurrence relation, explicitly given in Ref. 23. Introduc- 
ing a second function ry{u^s;q) {v=x, y or z), required 
while implementing the boundary conditions and defined 
in terms of Hy{u, s; q) according to 

d 

ry{u, s;q) = 1- fy{s; q){l -u^ 2q— In [Hy{u, s; q)]}, 

(2) 

where fx{s;q) = 1, fy{s;q) = I - sq and fz{s;q) = 1 - 
q, the polarizability of an ellipsoid with an anisotropic 
coating is 



Vo h{K=-l) 

47rng/„(/t= 1/ng-l)' 



(3) 



where 



-pvKiu-) + ^e^^][rl{u^) - SiS^^]. (4) 

where Si = Sin/^m^ and we have above introduced the 
short-hand notation rl{u±) = ry{u±,s; ef) and r^{u±) = 
ry{u±,s; el), is the standard depolarization factor for 
the 'u-direction for the outer surface and satisfies the sum 
^^^q23j24^. ^ _^^o haye also introduced 



Pv 



Ho{u+)Hl{u-y 



(5) 



where H^(u±) = Hy(u±,s',el) and Hl(u±) = 
Hy{u±^ s]el). We note that the total volume Vo only 
enters as a prefactor in the expression for a^y . The geom- 
etry of the particle enters through the geometric entities 
n^, rl{u±) and r'^{u±). The standard isotropic depolar- 
ization factor depends only on the shape {s and e^), 
whereas rl{u±) and r'^{u±) couple the geometry to the 
dielectric asymmetry of the coating. For an anisotropic 
coated sphere the expression for the polarizability agrees 
with the result obtained in Ref. 28. For an isotropic 
coating £\\ = the polarizability reduces to the stan- 
dard result given in for instance Ref. |^. The imaginary 
part of ayy is directly accessible through experimental 
extinction or absorption measurements. Explicitly, the 
extinction cross-section is Gyy = where ^[-J de- 

notes the imaginary part of the entity within the square 



brackets [] . In order to obtain the response of the coated 
metallic nanoparticle we must proceed by incorporating 
realistic microscopic dielectric functions for the coating 
and metal into the expression for the dipolar particle 
polarizability. 

Let us consider the interior metallic region. In a stan- 
dard fashion we assume that the dielectric function of the 
metal is described by a Drude function: 



(6) 



with Ljp being the plasmon frequency, Tp is a phenomeno- 
logical damping parameter and Soo is the dielectric func- 
tion for large frequencies. We take £00 = 1 throughout 
this study. An uncoated metallic ellipsoidal particle in 
vacuum (e^n = 1) has dipolar plasmon resonance frequen- 
cies 



(7) 



where nl is the (purely geometric) depolarization factor 
of the inner ellipsoidal surface^^*^^'^^ (see also appendix 

We now give the dielectric function of the coating; we 
assume that the coating consists of molecules with dipo- 
lar polarizabilities which are diagonal but with different 
components in normal and tangential directions; we de- 
note the polarizabilities by 7|| (7±), [where a subscript ± 
(II) denotes polarizability component perpendicular (par- 
allel) to the metallic surface normal]. Let us relate 7|| 
and 7^ to the dielectric functions Sc\\{oo) and €c±{^) ap- 
pearing in coated ellipsoid polarizability ayyi assuming 
that the surface is locally fiat, and imposing the con- 
ditions that the normal component of the total macro- 
scopic electric field and the tangential component of the 
displacements field are continuous across the surface sep- 
arating the molecular coating from the surrounding, we 
straightforwardly arrive at the following relation between 
the molecular polarizabilities and the coating dielectric 
functions2& 
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47r7x(cc;) 



^0 



e\\[u;) - 1 , 

Vo 



(8) 



where vq is the unit cell volume per molecule. Ex- 
plicitly, we take the following form for the renormalized 
molecular polarizabilities 



7m (^) = 7oo/i + (7o/i - 7oo/i)- 



(9) 



where 70^ is the static polarizability of the molecules, 
7oo/i is the high frequency polarizability, cjo/i is a res- 
onance frequency and is a damping parameter. 
It is sometimes convenient to characterize the electro- 
magnetic response properties of the coating by the ef- 
fective oscillator strength (compare to Ref. 16) cq^ = 
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47r(7o/2 — Too/ij/'^o, and the large frequency refractive in- 
dices n^ = ec±{oo = oo)"^^^ = (1 + 47r7oo/i/'^o)"^^^ and 

ny ec\\{u; 00)^^^ = (1 - 47r7oo/i/vo)""^^^. We point 
out that the induced dipole couphng between molecules 
(and image dipoles in the metal) in general, introduce 
extra anisotropy by renormalizing the resonance frequen- 
cies, static polarizabilities, large frequency polarizabili- 
ties and damping parameters compared to their "bare" 
values, ^^-^^-^^ and all such anisotropics are included in 
the expression 

To summarize, the general procedure for obtaining 
the polarizability avv{u) {v = x^y^z) for an ellipsoidal 
metallic nanoparticle coated with an anisotropic molecu- 
lar layer is: Consider a nanoparticle with principal semi- 
axes Qjx 5 CLy ■) and az. The size of the molecules together 
with these nanoparticle principal semiaxes determine the 
outer principal semiaxes bx, by and 6^, see Fig. Q] From 
these six principal semiaxes one calculates the the total 
volume Vo^ and the shape parameters s, e|, (for small 
coating thickness ef and are related through the rel- 
ative coating thickness Ay). The problem is completed 
by specifying the metallic nanoparticle parameters in Eq. 

and the renormalized polarizability parameters of 
Eq. ((2|); alternatively, one may use experimental results 
for the nanoparticle and coating dielectric functions. Fi- 
nally, the parameters above are used in the expression 
for ^^^(cc;) explicitly given in Eq. (0). 

For the case where the coating resonance frequency 
c^o/i is close to being degenerate with one of the parti- 
cle plasmon resonances (i.e., ujpy of Eq. Q) one expect 
hybridization between the two resonances.^ In order to 
address this point we proceed by finding an approximate 
expression for the resonance frequencies for ayy{uj). As- 
suming a thin coating Ay <C 1, and utilizing the near 
resonance approximation (the large i^-result^iSi) we ob- 
tain a simplified form for 0!,yfj , as detailed in appendix IXI 
Below we use this ayy for obtaining approximate expres- 
sions for the coupled (hybridized) resonance frequencies 
and for the extinction peak height at these resonances. 
The damping constants associated with the hybridized 
resonances are given in Eq. (|A28|) . We point out the the 
approximate expressions for ayy given in appendix IXI are 
useful for investigating other quantities as well. 

The resonance frequencies are obtained by finding the 
poles of ayy] using the results in the appendix we find 



that resonances occur at frequencies co - 
(see Eqs. (El2|) (IXT3|) . and (EEI)) 



]/2 and ujj^ 



]/2, with 



/ng; we have also introduced above a (dimen- 



sionless) parallel coupling strength 



nl/2 



Ay^^||Co||(l - <)- 



"oil 



(11) 



and a perpendicular coupling strength 



DyA 



Ay^^±Co±<- 



1/2 



(12) 



with the effective oscillator strength cqj^ = 47r(7o^ — 
7oo/i)/'^o as defined before, and the geometric factors de- 
fined as 



-'^Ojy 



and 



1 



-^ofy 



(13) 



1 



o2\-l 



where = fy{s;el) so that, = 1, 

= I- el and = 1 + (1 - sel)'^ + (1 - e^)" 
The quantities Qy^ appear also in the polarizability for 
an anisotropic shell^^ and satisfies the three sum rules 

9x\\^9y\\^9z\\ = gx±^gy±^gz± = 2 and^^x+^^ii = 1- 
Eq. ([TT]() has the same form as that for two dipole coupled 
oscillators, see for instance Ref. (also compare with 
Ref. |3). The strength of the molecular-plasmon interac- 
tion is characterized by the polarization and molecular 
orientation dependent coupling strength Dy^^ which de- 
pend on the relative coating volume parameter Ay, on 
the molecular response properties through co^i, as well as 
on different geometric quantities. 

For a sphere (e^ 0) we have that^y = 1/3, gy± = 
2/3, and = 1/3 and Eqs. CD and become 



D, 



v\\ Isph 



vl. |sph 




(14) 



The equations above quantify the degree of hybridization 
and its dependence on molecular orientation for a coated 
metallic sphered 

We notice from Eq. (PTTj) that if the coating and particle 
plasmon resonance frequencies are well separated, in the 
sense <C cj^^/cJq^, then the resonance frequencies 



are the uncoupled ones: cOpy and coq^ 



In contrast, for 
approximately at 



,4 



we have resonances 



V/l 



(10) 



D„ll « (l-0(Ay5„||Co||K)i/2, 



(15) 



(16) 



where for brevity, we have introduced the notations Thus for 



pv 



we get Dy 



and we have 



strong hybridization between the nanoparticle and coat- 
ing resonance frequencies, as will be illustrated in more 
detail in the next section. We also notice from Eq. ([T^ 
that for small n^, Dy\\ can be substantially larger than 
Dy±^ which can give a wider separation between u^^^ and 
compared with that for co^^. 



5 



The extinction peak height at one of the resonance 
frequency uof^ is (7yy{u = cj^^); using Eq. (|A7|) we 
convert this quantity into a scaled peak height = 
47r(7^^(a;^)/(Vo(a;^)^) at the resonances. Employing 
the thin coating approximation and assuming that we 
are close to a resonance Eqs. (|A14|) and (|A2Q|) apply, i.e. 
we have the approximate result 



UJZ, 



±\2 



which for 



pv 



■ uiQ^ gives 



1 



(17) 



(18) 



with Dy/^ given by Eq. (pj) . 

For the case of an uncoated metallic nanoparticle the 
scaled peak height takes the form Typ = l/{n^rp) (see 
Eq. (|A8|) ). Eq. (|TKj) shows that at strong coupling (i.e. 



for ujp^ ^ ^O/i), are of the same order of magnitude 
as Typ. 

For a very thin coating without any interior metal we 
have: Ty^{ei = 1) ^ AvgvfiCo^/T^ (see Eqs. (|Ali;) and 
(|A17|) . and Ref. 24), which is much smaller than Typ. 
Comparing with Eq. (|T7)) we notice that T^^/Ty^{£i = 
1) (X (Ay)~^ (note however that cj^^ depends on Ay 
), revealing substantial surface enhancement (see next 
section). 



III. RESULTS AND DISCUSSIONS. 

In the following, we illustrate the electromagnetic re- 
sponse properties of a coated metallic nanoparticle via 
the extinction cross-section ayy = uj^[ayy] for the case of 
a prolate {ay = a^, i.e. s = 1, "cigar-shaped") spheroid. 
We consider two different situations: the case when the 
molecules in the coating have their "resonant" axes (i) 
parallel and (ii) perpendicular to the surface normal. The 
polarizability component along the resonant axis is ay in 
case (i) and aj_ = a\\/2 in case (ii), which corresponds to 
coatings with identical but differently oriented molecules. 
For each of the two cases, the external field Eq is applied 
either along the long axis (x-axis) or the short axis {y- 
or z-axis). 

In Fig. the polarization averaged extinction cross- 
section for an uncoated metallic spheroid is shown to- 
gether with the cross-section for an spheroidally shaped 
molecular layer without an interior nanoparticle and with 
the resonant transition dipole moments oriented in the 
parallel and perpendicular directions. There are two 
plasmon peaks, as expected for a spheroidal metal par- 
ticle (two of the principal axes equal), with resonance 
wavelengths to the red (^o 1 1 ^- axis) and to the blue 



{EqWi/, z-axes) of the familiar LSP resonance of a sphere 
{X/Xp = a/3 ^ 1.73). The extinction peaks of the molec- 
ular layer is about 50 times weaker than the plasmons 
for both molecular orientations. Notice that, for geomet- 
ric reasons, the extinction peak height for the molecular 
layer is smallest for the case where the external field is in 
the x-direction, with the molecules oriented in the par- 
allel direction; see Ref. 24 for a thorough discussion on 
geometric and molecular orientation dependent effects in 
the response of an ellipsoidal molecular layer. Also no- 
tice the small difference in peak position between case (i) 
and (ii); in Ref. 24 it was shown that for Ay <C 1 there 
are no curvature induced shifts for the molecular layer 
resonance frequencies (see also appendix IX)) . Here, Ay 
is not sufficiently small and there are some minor shifts 
of the resonance frequencies. 

In Fig. 121 we show extinction spectra corresponding to 
the full coated spheroid response problem for the four 
distinct cases mentioned above, i.e molecular orienta- 
tion perpendicular (a, b) or parallel (c, d) to the sur- 
face normal and incident field parallel to the short (a, c) 
or the long (b, d) spheroid axes. The different colors of 
the curves corresponds to different coating resonance fre- 
quencies c^o/i- It is clear that there is strong hybridization 
between the coating and plasmon resonances whenever 
^Ofi ^ ^pv (m = 11?^ and V = x^y^z) in all cases. How- 
ever, the degree of hybridization differs greatly between 
the different molecular orientations and polarization con- 
figurations. 

In order to investigate the peak height of the reso- 
nances (see Fig. |2| Fig. |31show the scaled peak height 
(defined in Eq. (|A7|) ) compared to the scaled peak 
height Ty^{ei = 1) for a molecular layer (obtained from 
Fig. P), for different coating resonance frequencies. A 
comparison to the result given in Eq. ([T7)) is also made. 
We notice the two coupled resonances obtain similar peak 
height as c^o/i approaches the particle plasmon frequency 
ujpy. We also point out that the approximate expression 
Eq. (|T7)) works surprisingly well, considering that it is 
based on a close to resonance and thin coating approxi- 
mation, even far off the resonances. 

When the molecular resonance is far above or far be- 
low a plasmon resonance, we have a situation that can 
best be described as surface-enhanced absorption from 
the molecular layer i.e the magnitude of the molecu- 
lar resonance peak is greatly enhanced but its position 
and width is not changed dramatically from the case of a 
"free" molecular layer, see Figs |21 and [31 In this regime, 
the particle plasmons (at Ap^, with v = x^y^z) become 
either red-shifted or blue-shifted relative to the uncoated 
metallic spheroid resonance wavelength depending on the 
molecular resonance wavelength Aq . Interestingly, the en- 
hanced absorption from the molecular layer is not sym- 
metric on the two sides of the plasmon. This can be 
seen when comparing Fig. 2 a) and c) (see also Fig. 
corresponding to polarization along the short axis of the 
spheroid. When the molecular resonance is far to the red 
of the plasmon, absorption enhancement is seen for the 
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FIG. 2: Extinction spectrum for a prolate coated spheroid 
with different orientations of the molecules on the surface. 
The spectrum is in units of the total volume Vo of the coated 
spheroid, a^y — ix)'<s[avv\^ {/J^ = ||,-L and v = x,z) and 
A = 27rc/u; is the wavelength of the external electromagnetic 
field, a) The external field is along the z-axis (short axis) of 
the spheroid. The molecules have their resonant axes perpen- 
dicular to the metallic particle normal; b) the external field 
along the x-axis (long axis), and the molecules are perpendic- 
ular to the layer normal; c) the field is along the z-axis, and 
the molecules have their resonant axis parallel to the layer 
normal; d) the field is along the x-axis, and the molecules 
are oriented parallel to the layer normal. The different col- 
ors are spectrum for different molecular resonant wavelengths; 
the upper vertical bars indicate these resonance wavelengths 
Ao/x = 27vc/ooo^. Here, we chose the following values for ooq^: 
0.35a;p, 0.40;^, O.Gcjp, O.Tcjp, O.dup. The lower vertical bar is 
the particle plasmon resonance wavelengths Xpv = 27vc/cOpv 
{v = x,z). The same parameters as in Fig. Q]were used. No- 
tice the strong hybridization when Xq^ ^ Xpv and increased 
peak height (surface enhanced absorption) for the molecular 
resonance, compared with Fig. 



case when the molecular resonance axis is parallel to the 
surface normal (Fig. 2 c), but not for the perpendicular 
case (Fig. 2 a). The situation is reversed (although less 
pronounced) when the molecular resonance is to the blue 
of the plasmon, i.e absorption enhancement occurs for 
the perpendicular but not for the parallel case. Similar 
effects are seen in Fig. 2 b) and d), corresponding to po- 
larization along the long axis. From Fig. [3^) we notice 
that particularly strong surface enhancement (roughly a 
factor 50 for the extinction cross section for the param- 




FIG. 3: Scaled peak height for the coupled resonances in a 
coated prolate metal spheroid for the case a) when the coat- 
ing molecules have their resonant axis parallel to the surface 
normal and b) for a perpendicular resonance. The resonance 
frequencies oo^^ of the extinction (see Figs. |21and2J for difi'er- 
ent cjo/x was numerically determined and used to calculate the 
scaled peak height = 4navv{ootfj,) /{yo{(^tfi)'^) (y — x^z 
and /i = II, ^), see Eq. (|A7|) . The peak heights were normal- 
ized by the molecular layer peak heights T^^(£i = 1) obtained 
from Fig. Q The vertical solid lines correspond to the two 
plasmon resonances, ijjpx and cjp^, Eq. Q. The dashdotted 
horizontal lines correspond to the plasmon peak heights as 
given in the text right after Eq. (fT7)) . The solid and dashed 
curves are the approximate results given in Eq. (|17f) . The 
same parameters as in Fig. Q]were used. Notice the different 
ranges of the ordinates, and that the abscissas are in frequency 
units. 



eters used here) occur to the red of the plasmon for the 
parallel case with the incident field along the x-axis. 

The differences above can be understood as follows: To 
the red of a plasmon resonance, the induced field from 
the particle is in phase with the applied field and nor- 
mal to the surface at the poles of the particle (where the 
poles are defined by the direction of the induced dipole), 
resulting in enhanced absorption from molecules oriented 
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parallel the surface normal. Molecules with the perpen- 
dicular orientation instead couple mainly to the total field 
around the equator of the particle, where the field is per- 
pendicular to the surface normal. But to the red this 
field is weak, because here the induced field is out-of- 
phase with the applied field. Hence, there is very little 
enhanced absorption to the red of the plasmon in Fig 2 
a) but a large enhancement in Fig. 2 c). The same type 
of arguments apply to the case when the molecular reso- 
nance is to the blue of the plasmon, but here the situation 
becomes reversed because the induced particle dipole is 
out-of-phase with the incident field^ 

The case when the molecular resonance is degenerate 
with one of the plasmons corresponds to a regime that 
can not be described in terms of enhanced absorption. 
Instead, we have two completely hybridized resonances 
that exhibit " avoided crossing" , analogous to the case of 
two strongly coupled (quantum or classical) harmonic os- 
cillators, as is seen in Eq. Fig. ^illustrates this be- 
havior for the two molecular orientations considered. In 
the vocabulary of Ref. 3, the high frequency modes co^^ 
in Eq. correspond to "anti-bonding" combinations, 
i.e. the case where the induced dipole on the particle is 
out-of-phase with the molecular transition dipoles, while 
the "bonding" modes uj~^ corresponds to the case when 
the excitations are in phase. In a diabatic representa- 
tion, this corresponds to Rabi oscillations in which the 
excitation energy oscillates back and forth between the 
plasmonic particle and the molecular layer. As can be 
seen in Fig. ^ the mode splitting is much more signif- 
icant when the molecular resonance is oriented parallel 
to the surface normal (Fig. ^) than for the perpendic- 
ular case (Fig. This difference simply refiects the 
predominant polarization of the local field at the surface 
(the field would be strictly parallel to the surface normal 
for a perfect conductor at zero frequency). We also note 
that the splitting (expressed in frequency units) is larger 
for the plasmon that corresponds to polarization along 
the short axis of the prolate spheroid for both molecular 
orientations, which is somewhat surprising considering 
that the field enhancement is highest at the sharp ends 
of the spheroid^ However, the coupling strength is de- 
termined by the surface integrated local field at the two 
plasmon resonance frequencies, and for geometric reasons 
this quantity is higher for the doubly-degenerate short- 
axis plasmon in the present case. 

To the best of our knowledge, there are no experimen- 
tal studies that have directly probed the orientation de- 
pendent plasmon- molecule coupling discussed here. How- 
ever, the advanced nanofabrication and molecular func- 
tionalization technologies of today clearly make such 
studies a realistic possibility. One of the most interest- 
ing options could be to use a plasmonic nanoparticle cov- 
ered by an ordered layer of chromophores as a nanoscopic 
bio/chemo sensor. For example, it can be expected that 
the exact orientation of the chromophores on the sur- 
face will be sensitive to pH and to the interaction with 
molecules in solution. The chromophores might also be 
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FIG. 4: Illustration of strong coupling and "avoided crossing" 
between the plasmonic and molecular resonances of a coated 
prolate metal spheroid for the case a) when the molecules 
have their resonant axis parallel to the surface normal and 
b) for a perpendicular resonance. The horizontal dashed lines 
correspond to the two plasmon resonances, oopx and oopz, and 
the diagonal dashed lines denotes u = cjq/x- The solid white 
curves are the approximate results for the coupled resonances 
as given in Eq. The same parameters as in Fig. Qwere 

used. Notice that the axes are in frequency units. 

incorporated as functional groups into larger biomacro- 
molecules that change orientation/conformation upon 
biorecognition reaction with an analyte. Such sensing 
reactions would then affect the coupling between the sur- 
face plasmons and the chromophores, giving rice to pro- 
nounced changes in the extinction or scattering spectrum 
of the composite nanoparticle. To illustrate this possibil- 
ity. Fig. ISl shows the difference spectra between the par- 
allel and perpendicular orientations discussed above. As 
expected, the difference is largest for molecular resonance 
wavelengths that overlap the LSP modes. However, a 
noticeable difference is observed even far off the plasmon 
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FIG. 5: Difference in averaged extinction cross- section 
/\a /Vo — (cr" — (j^)/Vo (indicated by the different colors) 
between the parallel and perpendicular cases as shown in 
Fig. 2]a) and b), as a function molecular resonance frequency 
ujQ — — (jjQi. and external field frequency cj. Note that 
the axes are in frequency units. 



resonances, indicating that molecular orientation effects 
needs to be taken into account also in classical nanopar- 
ticle refractive index sensing experiments. 



APPENDIX A: APPROXIMATE EXPRESSION 
FOR COUPLED RESONANCES AND PEAK 
HEIGHTS AND WIDTHS 

In this appendix we use the large u (close to a reso- 
nance) and thin coating approximations for the coated 
ellipsoid polarizability given in Eq. in order to de- 
rive approximate expressions for the coupled resonance 
frequencies as well as for the peak heights and widths of 
these excitations. 

Close to a resonance {u large) and for a thin coating 
(Ay <C 1) we have that Eq. becomes 



Vo 



Jy{n -1) 



47rng Jy{n = 1/ng - 1)' 



with 



(Al) 



(A2) 



in an identical fashion as the derivation given in Ref. |2J. 
The geometric factors gy\\ and gv± {v = z) are given 
in Eq. (|Tt^ . Explicit expressions for the dielectric func- 
tions Si^ e± and e\\ appearing in Eq. (|A2|) are given in 
Eqs. Q and ((Hj). Below we investigate Jy(/^) in the limits 
of (1.) no coating Ay = 0; (2.) external field frequency 
close to the parallel resonance frequency, uj ^ ujq\\; (3.) 
external field frequency close to the perpendicular reso- 
nance frequency, co ^ u;o±. In (4.) we derive expressions 
for the peak widths of the hybridized resonances. 



1. No coating, Ay = 



IV. SUMMARY 

We have described an analytic method for calculating 
the optical response of arbitrary ellipsoidal nanoparticles 
with anisotropic molecular coatings in the small particle 
limit and applied this theory to the case of a prolate metal 
spheroid covered with chromophores. The results show 
that the hybridization between the molecular resonance 
and the localized plasmon resonances of the nanoparticle 
is highly anisotropic. It is suggested that this sensitivity 
can be utilized for novel bio/chemo sensing applications 
that are based on molecular orientation rather than re- 
fractive index contrast. 
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For the case of no coating. Ay = 0, we have that 
Jvif^) = Si-\- Hi and therefore the polarizability, Eq. (|A1|) . 
can be written 



47rng Ky{K = 1/ng - 1)' 



(A3) 



where Ky equals J„ multiplied a by a factor independent 
of K, and explicitly, 



(A4) 



with {^l-]) denoting the real (imaginary) part of the 
entity inside []. The resonance criterion ?R:[Ky{f<i = — 
1)] = gives, together with Eq. (|A4|) . the uncoated 
metallic ellipsoid resonance frequencies: 

^L,^ = KJ^ = (A5) 

in agreement with Eq. Q as it should (note that when 
Ay = then = n^). At one of the resonances the 
extinction peak height cr^^(cJres,^) = c^res,^^[Q^^^(^res,^)] 
is: 



Vo uj^[Ky{n = -I)] 



47rng %Ky{i<i = 1/ng - 1)] ' 
Vo 

_ I ■ 1 

reSjV 1 



-p-1, .2 



(A6) 
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where we have used Eq. ()A4|) and uj^es^v is given in Eq. 
(|A5|) . It is convenient to define a scaled peak height 



^^^^ 



SO that, using; Eq. ((SSl), we have 



(A7) 



(A8) 



Thus, for an uncoated nanoparticle Ty is a measure of 
the life time (inverse damping) of the plasmon excitation. 
We will below use the definition, Eq. (|A7|) . to compute 
approximate scaled peak heights of the coated metallic 
ellipsoid excitations. 



2. UJ ! 



' ^0|| 



Turning back to the general expression Eq. (rno|) we 
find that the resonance condition ^[Ky{i<i = l/n^ — 1)] = 
gives resonance frequencies: 

Kii) - ^ 



± 



1/2 



(A12) 



with a dimensionless coupling strength DyU defined by 



^vW = ^v9v\\Cq\\{1 - nl)[^ 

^^o|| 



+ ^Q|| 

2^011 



(A13) 



Using Eq. (|A7|) we find that the scaled peak height 
T^l (^^|) of the resonance characterized by resonance 

frequency uj\ (o;^, ) becomes (for ujp ^ 0) 



For the case uo ^ ujq\\ (and ujq \\ t^o±) we neglect 
the "perpendicular" term in Eq. \h2\ . so that Jy^K) = 
Ei -\- K,-\- AyEigy^e'^'^K,. Using the explicit expressions for 
the dielectric functions described by Eqs. 0, and (|S|) 
along with Q we find 



+ (cj^ + iu;Tp){u;'^ - cJq|| + iu;T\\)f^, (A9) 

Assuming that Ay <C 1, and neglecting terms of order 
(r = Fp, ry), and FAy (the neglected term Avgv±£± 
is of the same order as the terms neglected here) we find 
that the polarizability can be written in the same form 
as m Eq. (El now with 



-Av5„||Co||a;o||(Wp - uj^)k, 
+r^^{{u,l-u;\l + K))). (AlO) 



We notice that the second term in the expression for 
^[Ky{hz)] cannot be discarded because the first term may 
be close to zero. In the limit Ay,F|| the polariz- 
ability reduces to the uncoated result from the previous 
subsection. 

In the limit ujp^Tp (i.e. Si = Soo =13 recover. 



- p 

to lowest order in Ay, the result in Ref. 



q;„„(w) = Ay — 5.||Co|| ^ 



^11 



oil 



(All) 



T± = 



Ki-Kl)^) 



■r.K|-Kin+r||(c.2„g-(^±)2)' 



(A14) 



where we have used the fact that (^^|)^ is of the order 

(Ay)^/^, and neglected terms of order Ay. 

We notice from Eq. ^A12^ that if Dy\\ < \[ooln^ - 
a;Q||]/[2a;Q||] I, then the resonance frequencies are the un- 
coupled ones: cOp^/^ and cjoy- 

If instead D^ll > \ [u;ln^-u;^^^]/[2ujl^^]\, (i.e., forcj^n^ ^ 
^o||)' we get (cj^|)2 ^ ujln^[l±Dy\\], and 

(A15) 



-II ng(Fp + F||) 1 ±1^,11 

with Dy\\ = (Ay^^||Co||(l - O^/O^/^, showing strong 
hybridization of the plasmon and coating resonances. 



as it should. 



For the case co ~ cjqx (and cjoy 7^ ^o±) we neglect 
the "parallel" term in Eq. (|A2|) . so that Jy{K.) = Si -\- 
+ Ay^^^^x- Using the explicit expressions for the 
dielectric functions in Eqs. and along with Q an 
identical analysis as in the preceding subsection shows 
that polarizability can be written in the same form as in 
Eq. (El with 

^Ky{K)] = -u;(Tp{col^-u;^){l^i^) 

+F^(c^^-c^2(l + A.))). (A16) 
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In the limit Ay,r^ ^ we recover the uncoated result. 

In the limit oOp^Tp (i.e. = e^o = 1) we obtain, 
to lowest order in Ay, 



Vo 



47r uf 



■ — icc;r_L 



(A17) 



in agreement with the result in Ref. 124 

The resonance condition 5R[i^^(/^ = l/n^ — 1)] = gives 
resonance frequencies: 



± 



1/2 



(A18) 



with the dimensionless coupling strength Dy± defined by 

.2 „_...r->°+<y (A19) 



L>„_L = Avgv±co±n° 



compare to Eq. (|A13|) . 

The scaled peak height TJ^ (^Jl) of the resonance 
characterized by resonance frequency are ob- 

tained from Eq. We find (for cOp ^ 0) 



'0± 



in the same manner as in the preceding subsection. 



(A20) 



We notice now from Eq. (|A18|) that if Dy± <C 
]/[2a;Q^]|, then the resonance frequencies are the un- 



coupled ones: {jjp^Jri^ and a;o±. 
If instead D^^ ^ ll^p^v ~ 



2^o 



_), we get {uj^^Y 



: cc;^n^[l and 
1 



]/[2<]|, (i.e., for 



ng(rp + r^) i±i),j 



(A21) 



with Dy^ = (Ay^^xCox^^^)^^^, showing again strong hy- 
bridization of the plasmon and coating resonances. 



where Kv^\k,) = Ky{hi)\^^^± . Inserting this expansion 
into the expression for the polarizability Eq. (|A3|) . using 
the resonance condition 5R[i^^(/^ = l/n^ — 1)] =0, and 
keeping only terms to lowest order in Acj'^ we find 



.(0) 



?ft[Kir>{-l)]+i^Ki°\-l)] 



-. (A24) 



ALom[K',"il/nZ - 1)] + i9[ifr(l/ng - 1)] 

By furthermore neglecting the imaginary part in the nu- 
merator (i.e., assuming small damping) we find 



(A25) 



where we introduced the damping constant 



(A26) 



In order to obtain an explicit expression for the damping 
we proceed by expanding Eqs. lAlOjl and l|A16|l and 
obtain 

mi'\l/< - 1)] = ^U, - i^tf), (A27) 



where uj^^ = (cJq^ + cc;^n^)/2, and we neglected terms 
of order AyAcj^. Combining Eqs. (|A26|) and (|A27|) we 
finally find that the damping constant takes the form 



r±, - - + - ^ (A28) 



■ VjJ, 



>(<-(.;±)2) 



4. Damping constants 

In the preceding subsections we derived expression 
for the resonance frequencies and peak heights for the 
coating-nanoparticle hybridized resonances. We here 
proceed by also deriving expressions for the damping 
constants (inverse life times), which determine the peak 
widths. It is convenient to define 

Au,'^{u;t,f-u,\ (A22) 

wher e cj^^ are the resonance frequencies, see Eqs. (|A12|) 

and (IA18|) . We can then expand Ky{hi) in Eq. (IX3|) 
according to: 



(A23) 



Concluding, we above showed that close to one of the 
resonances the polarizability can be approximated by a 
Lorentzian, see Eqs. (|A25|) . with a width determined by 
Eq. (IX28ll . 

Using Eq. (IA24|) and the definition Eq. (IX7|) we find 
the following relation between the scaled peak height and 
the damping constant 



2(c.± )2 )2 r±^' 



(A29) 



to lowest order in Ay, compare Eq. I|A8(I . Inserting Eq. 
I|A28|) into the equation above we recover the results in 
Eqs. IIA14II and HA2nfl as it should. 
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